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^ . Abstract 

^ ' We calculate the neutrino mass matrix up to one loop order in the MSSM without imparity, 

! ■ including the bilinears in the mass insertion approximation. This introduces additional diagrams 

7— I ! usually neglected in the literature. We systematically consider the possible new diagrams, and 

find a few missing from our previous work. We provide analytic expressions for the mass matrix 
elements in the neutrino flavour basis, which are independent of the — Li basis choice in the 
I/"") , Lagrangian. We compare the contributions from different diagrams, and make "Lagrangian- 

basis-independent" estimates of when the new diagrams need to be included. We briefly discuss 
■ the phenomenology of a toy model of bilinear R-parity violation. 

t-H ! 

O ■ 

5 : 1 Introduction 

^ : 

T^" 1 ; The construction of a model which generates neutrino masses |IJ requires some extension of the 
^ Standard Model. Supersymmetry @, || is an interesting possibility which provides new (s)particles. 
y£^ • We consider the Minimal Supersymmetric Standard Model (MSSM) without imposing i?-parity. 
>- . The quantum number R p is denned as R p = (— \) L + 3B + 2S ^ where L, B, S are the lepton and baryon 
^ | number and the spin of the particle, respectively ||. In a general R p non-conserving supersymmetric 
5h I version of the SM, lepton number is not conserved^, which allows non-zero neutrino masses. 

The main physical motivation for neutrino masses comes from the anomalous neutrino data 
from both solar and atmospheric neutrino experiments ||. In the R-parity violating (]fl p ) MSSM, 
one neutrino can acquire a tree level mass through a see-saw effect from the mass matrix of the 
neutrinos and neutralinos ||. In order that the rest of the neutrinos become massive, loop diagrams 
that violate lepton number (by two units) must also be considered |7], || . For the model to generate 
neutrino masses that fit experimental constraints we must require that: Am^ ^ 1CT 4 eV 2 , Am 2 tm e 
[10 -3 , 10~ 2 ] eV 2 . For solar neutrino data there are several different possibilities for the relevant mixing 
angle (large or small mixing angle solutions) while for atmospheric neutrinos the mixing angle must 
be in the range sin 2 29 atm G [0.85, 1]. 

In the MSSM, the down-type Higgs and the lepton doublet superfields have the same gauge 
quantum numbers, which means that they can mix if lepton number is not conserved. Thus, we 
can construct a vector Lj = (Ha, Li) with J : 4..1. With this notation, the superpotential for the 
supersymmetric SM with R p violation can be written as 

\JKl 

W = fi J H u Lj + —LjL K E\ + X' J ^LjQ p D c g + h p t q H u Q p U c q . (1) 



1 Baryon number is also in general not conserved, however we will impose baryon number conservation by hand to 
avoid proton decay constraints, see e.g., |l 



We contract SU(2) doublets with e a p: en = 622 = 0, £12 = —£21 = —1- The R p violating and 
conserving coupling constants have been assembled into vectors and matrices in Lj space: we call 
the usual fi parameter /x 4 , and identify h J e k = A 4jfc ,Q and h p d q = A Apq . We call the usual ]ft p mass 
€i = Hi. Lower case roman indices i,j, k and p, q are lepton and quark generation indices. We 
frequently suppress the capitalised indices, writing Jl = (/i 4 , /i 3 , /z 2 , Hi)- We also include possible R p 
violating couplings among the soft SUSY breaking parameters, which can be written as 

V soft = ^HtH u + h^[ml] JK L K + B J H u Lj 

aJKI 

+ A ups H u Q p U c s + A Jps LjQ p D c s + —LjL K E^ + h.c. . (2) 

Note that we have absorbed the superpotential parameters into the A and B terms; e.g. we write 
B 4 H u H d not B 4 fi 4 H u H d . We abusively use capitals for superfields (as in eq. ([!])) and for their scalar 
components. 

In this notation where the Higgs joins the leptons as the fourth direction in {I//} space, the 
relative magnitude of R p violating and conserving coupling constants vary with the choice of the 
Higgs direction. R p violation can be understood geometrically: if all interactions in the Lagrangian 
agree which direction is the Higgs (e.g. H = L 4 ), then R p is a good symmetry. But if there is 
misalignment in {Lj} space between different couplings constants (for instance A ipq 7^ in the basis 
where Hi = /i 2 = ^3 = 0) then R p is not conserved. This Ifi p can be parametrised by basis independent 
combinations of coupling constants || |TD|, O, O. These invariant measures of R p violation in the 
Lagrangian are analogous to Jarlskog invariants which parametrise CP violation. 

Early work on lfL p neutrino masses f7| used models where the bilinear ]fi p could be rotated 
out of the mass terms (Bi, fii, rh^, refered to as "bilinears") into the trilinear interactions. Most 
subsequent analytic estimates have followed the calculation of and neglected the fi p bilinear 
masses in the loop contributions to the neutrino masses. This does not generate the complete set 
of one loop diagrams contributing to [mjy, so it is formally inconsistent. It would be acceptable if 
neglecting the bilinear contributions corresponded to neglecting loops suppressed by some additional 
small parameter. However, the size of the bilinears is basis dependent, so what is neglected depends 
on the basis choice, and this basis dependence of neutrino masses has caused some confusion in the 
literature. 

In a previous paper []I2"|| , we systematically analysed AL = 2 loop contributions to neutrino masses 



in the mass insertion approximation fHf . This led us to identify new diagrams which have not been 



included in many analytic estimates in the literature. We calculated the neutral loop diagram exactly 
(Grossman-Haber diagram), and gave basis-independent estimates for all diagrams using a common 
SUSY scale for all R-parity conserving mass parameters. The mass insertion approximation is valid 
as we know experimentally that neutrino masses are small. It is a particularly transparent way to 
include the bilinears, because the ]fl p Feynman rules can be read off the Lagrangian, and all the ^t p 
appears perturbatively in the numerator ( "upstairs" ) of each contribution. The relative contributions 
of different diagrams are easy to see. The main point of our previous paper was that the issue of the 
basis choice for calculating neutrino masses is transparent if we keep all contributions. 

The bilinear Jji p contribution to loops, as well as at tree level, has been included in [FJ, |T7 



where neutrino masses are obtained, after the exact diagonalization of the mass matrices of all 
particles which propagate in the AL = 2 loops. The bilinears mix Standard Model and SUSY 
particles, generating mass matrices of large dimension, so these results are partially numerical. They 
propagate tree-level mass eigenstates in their loop diagrams. The discussion in ref. |T7[ includes 



gauge loops. Reference have approximate diagonalisation formulae, calculated in the seesaw 



approximation, which is equivalent to the mass insertion approximation used here; they have few 



2 We have changed convention with respect to factors of 2 in A from our previous paper: we now put A/2 in the 
superpotential, with A 4y = h l J . 
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diagrams and long mixing angle formulae, whereas we have more diagrams and only MSSM mixing 
angles. 

The purpose of this paper is to provide basis-independent neutrino mass matrix elements, they 
are given in Appendices A and B. We present complete analytic expressions for all of the diagrams, 
and clarify how the neutrino mass matrix elements can be "basis-independent". We discuss in 
detail when the different contributions to the neutrino mass matrix are relevant, identifying for both 
bilinear and trilinear contributions when they must be included or can be neglected. We discuss 
phenomenological issues for a toy model where the R p violation comes from the soft SUSY breaking 
sector of the potential. We will discuss phenomenology in more detail in a subsequent publication 
7M. In the literature, most analyses using neutrino data have taken into account only the effects of 



Rp violation from trilinears and the fj, bilinear [18 . 

In section |2], we discuss what we mean by "basis independent" neutrino mass matrix elements, 
what we calculate, and how basis dependent it is. In section [3], we discuss when the bilinears 
are important in the loop contributions to neutrino masses. This is a detailed discussion, with 
generation indices, of results presented in \T2\. In section [|, we study the phenomenology of some 



simple bilinear models. We conclude in section ^. Exact one-loop neutrino mass matrix elements are 
presented in Appendix A in the pi) = basis. In Appendix B, the neutrino mass matrix elements 
are given in terms of MSSM parameters and basis-independent combinations of coupling constants 
that parametrise Jft p . Additional Feynman rules for including Jft p in the mass insertion approximation 
can be found in Appendix [D]. Numerical bounds on R-parity violating parameters obtained from the 
different diagrams are given in Appendix 



2 Issues of basis 

There is no unique interaction eigenstate basis for the Y = — 1 fields (the components of the Lj 
and Hd superfields) in the Tji v MSSM. This means that the Lagrangian parameters depend on an 
arbitrary choice of basis, or equivalently that Lagrangians which differ by a rotation in \L{\ space 
make equivalent physical predictions. Observables cannot depend on the choice of basis in the 
Lagrangian C, so they must be scalar functions of the vector and tensor parameters in {Lj} space. 
It is common in SUSY to study the dependence of physical observables on Lagrangian inputs. The 
basis-dependence of the Lagrangian makes this problematic in ]fl p models because a point in physical 
parameter space corresponds to different numerical inputs in different bases. So comparing results 
calculated in different bases is difficult. 

The lfl p in the Lagrangian can be parametrised in a basis independent way by combining coupling 
constants, masses and vacuum expectation values (vectors and matrices in {Lj} space) into scalar 
"invariants" , which are zero when R p is conserved. Physical observables can be expressed in terms 
of these invariants and other R p conserving quantities. For a discussion of constructing invariants, 
see e.g. Ml O O, 0. 



Neutrino masses and mixing angles are measurable quantities, so must be (Lagrangian) basis- 
independent. The neutrino mass matrix [m„]jj, i,j : 1..3, is basis dependent, depending on the 
choice of Higgs direction in {Lj} space, and on the choice of lepton flavour directions Lj and Lj. 
A physically motivated basis choice for the Lj is the charged lepton mass eigenstate basis (= the 
neutrino flavour basis), because the neutrino masses and mixing angles can be computed from [m^jj 
in this basis, so it contains all the measurable information. We can express the charged lepton mass 
eigenstate vectors in a basis-independent way, using the mass insertion approximation, so we present 
Lagrangian basis independent expressions for the elements of [m v ]ij in the flavour basis. So we have 
"basis-independent" expressions for matrix elements which are intrinsically basis-dependent — this is 
sensible because we want expressions which are independent of the arbitrary Lagrangian basis choice; 
the flavour basis is "physical" . 
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Table 1: The basis-independent invariants used to parametrise the ftp relevant for neutrino masses, 
together with their value in the (p) = basis. They are zero if R p is conserved. 5^ and 5b parametrise 
bilinear Jjt v . Note that these invariants have signs: for arbitrary vectors a and b, a ■ X 1 ■ b = —b • A 1 • a. 



The basis we work in is 



V 



H d = - 
|A i • v\ 



(3) 



where v is the vector of vacuum expectation values in the down-type Higgs sector and we impose 
the requirement that v ■ \ l ■ V ■ v oc 5^. The lower case index labelling the matrix X k is the 
singlet lepton index: [\ k ]u = \ IJk . So there are three flavour eigenstates L l in {Lj} space: Lj = 
\ JKi vk / \/vp\ NPi \ NMt VM (no sum on i) 0. This corresponds to the charged lepton mass eigenstate 
basis in the absence of In this R„ conserving case A 4 *- 5 = h l J and \ kl i = 0, where k, i,j : 1..3. So 



\ Uk Vl 



\[2m? k (I 



4, J 



X 1 ■ X 1 ■ v oc 5 V is the condition that the {e^.} are in the mass 



eigenstate basis, and A* ■ v are the charged doublet eigenvectors. In the presence of bilinear ftp , the 
basis (H) will not be exactly the charged lepton mass eigenstate basis. However, bilinear ^ p masses 
are required to be small by neutrino masses, so the basis (|3|) is close to the charged lepton mass 
eigenstate basis. We include the correction due to the bilinears in the mass insertion approximation. 

Neutrino masses are proportional to Jfip couplings, and observed to be small. We can therefore 
compute neutrino masses by perturbing in Jjl v masses and trilinear couplings. We include the 
bilinear Jji v masses in the mass insertion approximation, which means all lft p couplings appear in the 
numerator ("upstairs") of the expression corresponding to a diagram. We propagate MSSM mass 
eigenstates, whose masses we make "basis-independent" by substituting _B 4 — > B-v/\v\, \i± — > j2-v/\v\, 
and rn.44 ~^ v ' ^ ' ^/|"^| 2 - The physically relevant flavour basis is written in terms of vevs and 
coupling constants (see equation |3|), so neutrino mass matrix elements in this basis are proportional 
to invariants. We list the four invariants relevant for neutrino masses in table 0, along with their 
value in the basis where the sneutrinos have no vev. 

A significant advantage of the basis-independent formalism is that the neutrino masses can be 
computed in any Lagrangian basis, including those where the mass insertion approximation is not 
valid. This could be easier than rotating to the (z/j) = basis, because the vev is a derived quantity, 
calculated by minimising the potential. For instance, consider a model where all the Jft p is in B and 



3 The capitalised index ordering on A ' is because o ■ 
of L l is -v- X l /\v- \% 



A ■ b = —b ■ A • a is an antisymmetric product. The transpose 



A 



v — that is, there is a basis where all the Jft p couplings other than Bi and (pi) are zero. Clearly it is 
easier to evaluate invariants in this basis than to rotate to (Pi) = 0. The basis independent formalism 
also makes it simple to compare results computed in different bases. 

The tree level neutrino mass is non-zero if 5^ ^ 0. The diagram appears in figure [3] in the 
mass insertion approximation. In the basis (|3|), it contributes a mass matrix 

a m Xa 

which gives a mass mf 66 = £i, a (^u) 2 Z^\[i\ 2 /rn Xa to the neutrino 

vT e = Y-Li, (5) 



where 5^ = yX^^) 2 ; an d the neutralino index a runs from 1 to 4. The index "3" on Z corresponds 

to the interaction eigenstate h u — see Appendix C for our conventions on Z. Equation (^) is equivalent 
to the usual formula m*J' ee = det[M^]/det[M^}, where (JlfW) is the % 5 x 5 neutralino and 
neutrino mass matrix (MSSM neutralino mass matrix), as can be seen by writing and in 
the MSSM mass eigenstate basis with (Pi) = 0. 
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Figure 1: Tree- level neutrino mass in the mass insertion approximation 

Loop corrections to the neutrino mass matrix can be divided into three categories. First, there 
are gauge and top Yukawa coupling loops which renormalise the mass of z/g ree — we neglect these 
because their effect is small |15| . Secondly, there are loop corrections to 5 l (equivalently, \ii and/or 



(Pi)) which can modify the direction of P3 — these we partially compute and list in the Appendices. 
They are renormalisation scale dependent, because they are loop corrections to the tree mass. We 
are not interested in loop corrections to the tree mass, so we rotate these away. This is discussed 
in more detail in Appendix |A]. Finally, there are finite loops which give mass to any neutrino. The 
third group are the most interesting loops, because they generate mass for the two neutrinos who 
are massless at tree level. Schematic representations of the one-loop neutrino mass diagrams that 
we consider are reproduced in figure [|. This is a slightly modified version of a figure from ref. [|I~2jl . 
Each of the four diagrams represents a number of Feynman diagrams. This should be an almost^ 
complete set of one loop, AL = 2 diagrams which generate mass matrix elements for the neutrinos 
that are massless at tree level. The possible places on the diagram for the two required lepton number 
violating interactions are labelled I.. VII I. With the definition H oc v , lepton number violation is 
not possible at the charged lepton mass insertion in diagrams a) and d). (A sneutrino vev could 
provide lepton number violation at this point in a basis where (Pi) ^ 0.) 

In Appendix ||, we list the loop diagrams, giving exact formulae in the (Pi) = basis for each 
diagram's contribution to [m v ]ij. The basis-independent version of these formulae is in Appendix 
y. In table 0, we summarize the diagrams and make basis independent estimates by setting all the 
heavy masses to a unique scale rrtsusY- The lists starts with the canonical trilinear diagrams, then 
the Grossman-Haber loop induced by the soft bilinears, and finally all the additional diagrams which 



We discuss in Appendix ^| why we neglect certain finite diagrams of the third type. 
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Table 2: Estimated contributions to [m^ from all the diagrams. In the second two columns is the 
label of the diagram of figure 2, and the position on the diagram of the two AL = 1 interactions. 
Column four is the "basis independent" estimated contribution to the neutrino mass matrix in the 
flavour basis. All indices other than i and j are summed. 



arise when the bilinear 8^ are included in the loops. For each diagram, labelled by a — > d for the 
diagram category of figure 0, and by roman numeral identifying where the ^t p should appear on the 
figure, we give an estimate of the diagram in terms of the invariants listed in table |l|. We have added 
a few relevant diagrams missing from the original list. In Appendix B, we make an exhaustive list 
of all potential diagrams, explain why some are zero, and give basis-independent expressions. 

To evaluate diagram b) with Iji v at points // and VII (this is the usual coloured trilinear diagram), 
we identify the external neutrino legs as Lj and Lj. At vertices II and VII sit couplings A sp and 
\' rs . (The bold face indices correspond to the squark mass eigenstate basis — see Appendix [D].) The 
diagram is therefore 



m 



u\kj 
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s,p,r 

35 k x f p 5{r \m ds m dp \ 
8tt 2 msusr 

(8) 

where A? r = -[(A\') Rpr v R + v u fi R \' Rpr ], and we work in the down quark mass eigenstate basis: 
X' st ■ v is diagonal. The integrals / are listed in Appendix |B|. The last line corresponds to the last 
column of table |2|. Note that we do not divide by 2 when we symmetrise on i and j, because the 
diagram with % <-> j is different. This agrees with [[16], [17] f] and disagrees with ||. It is clear in the 
mass insertion approximation that one should not divide by 2, because in one case V{ couples to D c 
and in the other case to Q. In the mass eigenstate formalism, one can see that these are distinct 
diagrams not included in the sum over p, r, s by considering the case where only A 132 and A 223 are 
non-zero. 

In our previous paper, we were unclear about whether the charged goldstone boson of SU(2) 
could mix with the {E c }. This was due to various sign discrepancies in the literature. As expected, 
the goldstone is pure SU(2) — so there is no loop diagram propagating a W 11 that is proportional to 
a gauge times a Yukawa coupling. 



3 When are the bilinears important in loops? 

Bilinear contributions to loops have traditionally been neglected, which intuitively seems reasonable 
if they are "small" in a "sensible" basis close to the MSSM. However, this apparently reasonable 
oversight is confusing, because the size of what is neglected depends on the basis choice. In this 
section, we provide basis-independent conditions for when the bilinears can be neglected in the 
loops — assuming the mass matrix elements are then calculated in a sensible basis where the bilinears 
are small. These are not hard and fast rules, but estimates of when the extra bilinear diagrams should 
be included. 

A slightly different approach would be to consider the relative size of all the diagrams, and 
catalogue all the permutations of which diagrams should be included and which neglected. For 
instance, it is possible that the trilinears are negligeable with respect to the bilinear diagrams, or 
that the A trilinear contributes only a small correction to the A' trilinear, and can be ignored. Here 
we assume that the trilinear contributions are always calculated, because they take little effort. Then 
we ask whether 5 b should also be included — this is a few more diagrams. Finally we consider adding 
the 5^ contributions, which is many more diagrams. We do not specifically discuss the case where 5^ 
should be included but the 5b contributions are insignificant, because the 5b diagrams are relatively 
little work compared to the many 5^ diagrams. 

In our previous paper, we outlined three cases: 

• A: all the bilinear contributions to loop diagrams are negligeable. In this case the loop contri- 
butions to the neutrino mass matrix are the usual trilinear diagrams of figures 2a and 2b. 

• B: 5^ is negligeably small, but 5b should be included. In addition to the diagrams of case A, 
one should consider the neutral Grossman-Haber loop of figure 2c, and possibly the 5b mass 
insertion to figure 2a. 

• C: Include all bilinear JfL p contributions. There are additional contributions to diagrams 2a, 
2b, and 2c, and a new diagram 2d. 

5 We thank E.J. Chun for a discussion of this point. 
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We determine basis-independent criteria for when the bilinears can be neglected in the loops by 
comparing our estimates of the size of each diagram from the last column of table §. A diagram 
should be included if its contribution is of order the trilinear loops — however, we want to avoid 
comparing bilinear loop corrections of the tree mass (which are irrelevant), with the trilinear loops 
contributing masses to neutrinos which are massless at tree level. So we distinguish two cases : 

m tree < m loop ^ and m tree > m loop _ 

If m tree ^ m loop , we simply compare bilinear to trilinear loop mass matrix elements. Since 5^ is 
not large enough to produce m tree > m l ij° p , it is negligeable in the loops, so we are in cases A or B. 

If m tree 3> m loop , we consider only the loop mass matrix for the neutrinos that are massless at 
tree level. We compare the bilinear and trilinear contributions to this sub-matrix. Cases A, B, or C 
can arise if w^ ee 3> m loop . 



3.1 case B — including the soft bilinear 

We first consider when 5b should be included in the loops. The Grossman-Haber loop is of order[] 



g 2 5 t B 5 B m x /(64:TT 2 ), which is potentially large because it is proportional to gauge rather than 



Yukawa couplings. Recall that the canonical trilinear diagrams are of order m l J ~ 5y 3 5^ 3 ml/(87r 2 msusY)- 

If the tree mass is small, m tree ^ m loop , then the Grossman-Haber loop should be included if it is 
of order the trilinear loop. This will occur if 

9 -si > s{fy/m, 5{ ki jm, (9) 

or the condition to neglect 5b can be roughly estimated as 

5 j B < 5$ 3 h b , 5{ 33 h T (m tree < m loop ), (10) 

Alternatively, if m tree 3> m loop , we are only interested in the loop contributions to the mass of 
neutrinos z/ 2 and u% who are massless at tree level: if B is aligned with ft, then the GH loop is 
a correction to the tree-level mass, and can be neglected along with other such loops. B can be 
decomposed in components parallel and perpendicular to //: 

B = B ± + B h (11) 

with B\\ = (B ■ jl)j2/\j2\ 2 . The part of the Grossman-Haber loop proportional to 5 B ^5 B ^ is a loop 

correction to the tree mass m tree . The 5 B ^5 B± terms mix v tree with v 2 and v\ y and the mass matrix 
from the GH loop for v\ and v 2 is 

This will be comparable or greater than the trilinear loops when 

fo ± = f & - -^k) a sf^/m- (is) 

— * op 

We define 5y P = 5y P — S^JX' ■ p)/\fl\. This translates roughly into the condition that 5b can be 
neglected in the loops if 

5 Bi <^5§ 3 h b ,5f 3 h T K ree »m^). (14) 



6 We are neglecting flavour violation among the sneutrinos, which could induce a significant loop mass even if 
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3.2 case C — including [i{ 



The tree-level neutrino mass is ~ Z)i(^u) 2 ' 7l x> so ^ $1 * s re l evan t m the loops, then m loop <C m tree . 
The additional diagrams proportional to <5^ should be included if they induce a loop mass for v 2 and 
z/i greater than or of order the trilinear or GH loops. 

We first consider the various 5^ contributions to diagrams a and d of figure |2] which are of order 
~ [^u(h % e ) 2 ][Si(hP e ) 2 ]msusY- Generically the vector (5j 1 h T h T , d^h^h^, 5pi e h e ) will not be aligned with 
(8*, 6^, 5*), so the 5^ bilinear contributions should be included if 



^(K) 2 ~ S x V h eKi V h d h % h j, k, I, p, q not summed. (15) 
A rough guide to when these 5^ corrections can be neglected is therefore 

#/r«ar,aW£). (i6) 

There are also bilinear loop contributions of the form 

W oop h~ %{KfWMsusY (i,jG{r jA *,e}). (17) 

These initially appear more significant, because they are only suppressed by two, rather than four, 
trilinears or yukawas. However, m l ° op <C mjj" ee and v tree oc (5*, 8fc, 5^), so contributions of the form of 
equation ( |TTD mix tree with v 2 and v 1 . These latter two neutrinos, which are massless at tree level, 
therefore acquire a seesaw mass of order 



k 



~ K(Kr}[(K)\}rnsusY (i,je {3,2,1}). (18) 

This is the mass matrix structure we considered before equation ([15|) , and will make a significant 
contribution if equation ( jT5|) is satisfied. So equation fll5"D is the condition for when 5^ bilinears 
should be included in loops, or equivalently, equation ( |16"D is a rough estimate of when the 5^ loops 
can be neglected. 



4 Phenomenology 

In this section we briefly consider the constraints that can be set on Jji v couplings from solar and 



atmospheric neutrino data. We follow the approach of p0| , 21| , who set bounds on ]fl p models 
corresponding to case A — models where the bilinears are negligeable in loops and neutrino loop 
masses are due to trilinear couplings. We discuss here a toy model where the is in the bilinears, 
and the 5b loops make a significant contribution to \m v ]ij (case B of the previous section). We leave 
case C — where 5^ should be included in the loops — for a subsequent analysis plf . We would like to 
establish bounds on the different set of basis independent combinations of coupling constants that 
contribute to the neutrino mass matrix. These bounds are more indicative of orders of magnitude 
than steadfast constraints. In order to proceed we make the following assumptions. First, we use the 
results of reference which constrain a general 3x3 symmetric mass matrix from neutrino data 
considering: atmospheric and solar neutrino experiments, CHOOZ experiment and the constraint 
from neutrinoless double beta decay. The overall conclusion of this analysis is that the mass matrix 
elements should be constrained to be on the order of | f^]^! ^ O.leV. Here we present our bounds 
allowing each individual matrix element to take on the maximum value of 0.1 eV. This simple 
approach can be extended to obtain bounds when we include combinations of contributions to the 



o 



neutrino mass matrix from 5* 5 l B , 5 l £ ' , for any given model which contains these R-parity violating 
terms originally in the Lagrangian. 

Secondly, we assume that each of the contributions from the diagrams that we have found is 
separately the unique contribution at a given timeQ. In this way we are trying to obtain the weakest 
constraint on the basis-independent combination of coupling constants. That is we allow each sepa- 
rate contribution to take on its maximum value. We can then choose from all of the bounds which 
is the strongest constraint on a given combinations of coupling constants. We also assume, unless 
otherwise indicated, that all supersymmetric mass scales are the same of order m snsy . The next 
approximation we make is that when we sum over neutralinos or charginos there is no suppression 
arising from the mixing angles. We also assume that the magnitude of all R-parity violating couplings 
constants to be generation blind separately. That is, all Bi to be of the same order of magnitude, 
and similarly for /ij, Ay&, A^- fc , but we do not suppose beforehand that R-parity violating coupling 
constants arising from different terms in the Lagrangian are also of the same order of magnitude. 
The numerical results for all diagrams are summarized in appendix [E[ 

As we can see from the tables in appendix |E] the strongest bounds for the combination of basis- 
independent coupling constants S B S B arises from the Grossman-Haber diagram: 

5 B 5 J B < 3 x KT 10 . (19) 

For 5^5^ the best bound comes from diagram 19, which is the neutral Grossman-Haber loop with 
R-parity violation at points I and VI: 

4^ < 7 x KT 10 . (20) 

(We use tan/5 = 2 and a generic m susy ~ 100 GeV for both these bounds.) The strongest constraint 
on 5^5^ is from the tree-level mass: 

^<10" 12 • (21) 

For comparison, we list the bounds which can be derived from the remaining diagrams in Appendix 
??. 



4.1 Toy Model 

Suppose we only allow the presence of the basis-independent combinations of coupling constants 5^ 
and 5 B . A non-zero value of 5^ can arise either from having /ij ^ or < z>j >^ 0, so this model 
could arise if all the lji v originates in the soft SUSY breaking terms Bi, which induces a misalignment 
between the vev and superpotential couplings. (This would induce 5y ~ hbS^ and 5\ ~ h T 5^, which 
we neglect because we assume our model to be in case B, where contributions of this size can be 



neglected.) This differs slightly from the usual bilinear model, discussed in detail in [To], [L6[ |Ty| |, 
where the lji v originates in the GUT-scale misalignment between pi and the trilinears, so that B 
becomes misaligned with respect to fl while running down to the weak scale. From the results of 
appendix [E| we see that the relevant contributions to the neutrino mass matrix will arise from the 
tree level contribution plus diagrams 3 and 19. The upper bound 5^ ^ 10~ 6 from the tree level mass 
contribution ensures that (for low tan/3) the remaining diagrams, which in principle contribute to 
the neutrino mass matrix through mass insertions of 5^ in the loops, are negligeable. 

With these three contributions we can obtain neutrino masses in the phenomenologically interest- 
ing region which can simultaneously satisfy atmospheric and solar neutrino data constraints. There 
are several ways to see this. The full mass matrix using our simple approximations is given by, 

rn$ = m susy <^ + axm SWY 8 B 8 B + a 1 m susy (^^ + 5^), (22) 



7 This may be incorrect in a poorly chosen basis, where different diagrams cancel each other UA 
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where a% ~ gf - ?- Certain subcases of this mass matrix can generate neutrino mass textures, which 
have been analyzed in the literature, that can produce neutrino masses and mixing angles in accor- 
dance with present neutrino data. As a simple example, we see that if we set 5% = 0, and neglect 
the S^6b term in equation (p2|), then the structure of the neutrino mass matrix corresponds to the 
one analysed in case 1 of ref. pl| . In this case the neutrino mass matrix at one- loop consists simply 
of the tree-level term plus a correction to the 1-2 submatrix. The subcase analyzed in ref. |2T| had 
this same structure but the loop correction was due to the trilinear diagram. 

It is well known that the tree level contribution can give mass to only one of the neutrinos. 
Together with the loop contributions one of the neutrinos will have a mass 



77lj,3 ^susy 



El^ + v^iJ 2 , (23) 



where we have used the decomposition B = Bj_ + B» of equation (|TT|). 

The sub-mass-matrix for the remaining two neutrinos is 1 and is 2 is given by 

ml j = aim susy 4 x 4 x - (24) 

Terms of the form S B± 5^ mix the heavier is 3 and the lighter (is 1 , is 2 ). It can be seen from equation 
(|18|) that the seesaw mass generated from this mixing is suppressed, given our choice of SUSY 
parameters, compared to the direct loop contribution given in eq. ([24]). So in this case, this model 
corresponds to case 1 of pi ]. This simple example f\ shows that a hierarchical neutrino mass spectrum, 



which satisfies atmospheric and solar neutrino data constraints, can easily be obtained by taking 
|<5/i + v^i^Bn I ~ f ew x 10~ 6 and 5 l B± ~ few x 1CT 6 — 1CT 5 . This gives us a massless neutrino and two 
massive neutrinos such that Am ata ~ m 2 3 and Am so i ar ~ m 2 2 - It is also possible to obtain other 
types of spectra for different input values. 



5 Conclusions 

There are many diagrams which can contribute to the neutrino mass matrix in the framework of the 
MSSM without R-parity. In general, in the literature only a few of these have been considered. In the 
present paper we have obtained the full analytic expression for the different diagrams contributing to 
the neutrino mass matrix from R-parity violating bilinear and trilinear coupling constants. We have 
expressed each contribution both in the (vi) = basis and in terms basis- independent combinations 
of couplings. We have also shown when the separate diagrams should be included/are relevant in 
the context of a given consistent framework. We have presented bounds on combinations of the 
basis-independent couplings constants from neutrino experimental data, and shown that a simple 
toy model of bilinear R-parity violation can successfully accomodate the necessary mass squared 
differences to account for the neutrino oscillations. 



Acknowledgements 

We would like to thank Eung Jin Chun, Herbi Dreiner and Paolo Gambino for useful conversations. 
The work of M.L. was partially suported by Colciencias-BID, under contract no. 120-2000. 



Appendices 

In these Appendices, we present the one loop contributions to the neutrino mass matrix, excluding 
the one-loop corrections to the tree level mass, because we are interested in the one-loop masses for 

8 We leave a more detailed numerical analysis for future work. 
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the neutrinos that are massless at tree level. We identify these loops in Appendix [A]. In Appendix 
[B], we list diagrams in the order of the table 2. We give the amplitude for each non-zero diagram in 
the (pi) = basis. In Appendix we systematically list all diagrams, ordering them by the roman 
numerals that identify in figure |2| where are the two units of Jji v . The amplitudes in Appendix B 
are "basis-independent", that is they are expressed in terms of invariants and MSSM parameters. 
Notice that the neutrino mass matrix elements are minus the amplitude for the diagram, [m u ]ij = 
—M.ij{p 2 = 0), so our formulae are for — [mjjj. We present the ]fi p Feynman rules in Appendix |D[ 
Appendix [E] contains our numerical results for bounds placed on combinations of basis-independent 
couplings constants from neutrino data. 



A renormalisation — which are the finite loops 

We wish to neglect loop corrections to the tree level mass, because we do not need such accuracy and 
because we prefer to avoid the issue of renormalisation. We therefore neglect all diagrams involving 
gauge bosons and those with ]fi p at I and VIII. We initially expected the remaining loops to be finite, 
because they could contribute mass to the neutrinos who are massless at tree level. However, diagram 
17 (figure 2d with ft p at I and VII) is renormalisation scale (Q 2 ) dependent. In this Appendix, we 
show that the offending Q 2 dependent diagram is a loop correction to the tree mass. 
The tree mass matrix can be written 

K]y = mT\eT e WT% (25) 

where e| ree is the eigenvector associated with m l { ee . At one loop, both the mass and eigenvector are 
modified: 

K] fi = « ee + Am 3 )(e* ree + ^hWT + Ae 3 )i (26) 

We do not want to include loops contributing to Am 3 or Ae 3 . We do not calculate gauge loops and 
diagrams with Jji v at I and VIII, which contribute to Am 3 . However, there are also loop corrections 
which change the direction of e 3 , that is loop corrections to the angles in the rotation matrix which 
diagonalises the 7x7 neutral fermion mass matrix. These loops are included in our calculation, and 
we want to identify them and throw them out. 

We work in the mass insertion approximation, in the flavour basis where (Li) J oc A This 
is the basis where the charged lepton mass matrix is diagonal in the MSSM. Some care is required 
in determining the {vj}. For one-loop neutrino masses, it is sufficient to use the vevs {vj} that 
minimise the tree level potential. To identify the one-loop corrections to the tree mass, we must use 
some one-loop choice for the {v j}. There are various possibilities, for instance we could use the vevs 
which minimise the one-loop effective potential, or we could use the tree + one-loop masses that mix 
neutrinos and h° d with gauginos. We opt for the latter, because in such a basis it is easy to separate 
the loop corrections to m^e^e^ from the loop masses m 2 and m\. So we choose a basis where there 
are no tree or one- loop mass terms mixing w° with v^. This means that there will be small tree- level 
sneutrino vevs which cancel the one-loop w°Vi mass. These vevs are formally of one loop order, so are 
irrelevant inside the loops, because their contribution would be of two loop order. The contribution 
of these sneutrino vevs to the numerical value of the tree level mass is also a higher order effect and 
therefore negligeable. The one place these vevs must be included is in the one-loop z/jfe mass, where 
the tree level sneutrino vev contribution (formally one loop order) cancels the scale dependent part 
of the one-loop v^b mass, but leaves a finite v^b mixing. These can contribute to the loop masses 
m 2 and mi, as can loop contributions to h^Vi mixing. 

We know that in the flavour basis at tree level the neutrino mass is given by eqn [|. The one-loop 
expression will be the same, in our present basis, provided we identify the ^ of eqn [| with the 
tree+ one- loop mass terms that mix h° u with v^. Since we are only interested in the lowest order 
contributions to neutrino masses, we only include the tree level contribution to m 3 . Knowing the 



one-loop corrections allows us to identify the finite loops that generate masses m 2 and m\. these will 
be loops without mass insertions on the external legs, and loops with one mass insertion that mixes 
a neutrino with ha or b°. The loops mixing a neutrino with hd or b° are finite in our present basis, 
as we show in Appendix ||, diagram 17. 

In summary, we identify the loops contributing to one-loop neutrino masses m 2 and mi to be 
those without mass insertions on the external legs, and diagrams with a mass insertion on one of 
the legs that mixes a neutrino with the b° or h° d components of a neutralino. The remaining loop 
diagrams, with a mass insertion on both external legs, or a mass insertion mixing z/j with w° or h u 
are corrections to the tree mass. 



B Diagrams in (fit) = basis 

We define 



.4 



ml 



and 



.4 



ml 



« iml 

d 



(27) 



:28) 



We include A and Ad in the mass insertion approximation, and seperately diagonalise the 3x3 
slepton ( or down squark) doublet and singlet mass matrices. Bold face indices indicate these slepton 
or squark mass eigenstate bases, see Appendix C for details. Lepton and down quark mass eigenstate 
indices are in ordinary type. 

We also define 



J(m 1 ,m 2 ) 



1 



m; 



._ l n ^i. 

167r 2 m\ — m\ m\ 



(29) 



I(m 1 ,m 2 ,m 3 ) 



d 4 k 



(2vr) 4 k 2 + m\k 2 + m\ k 2 + m\ 



,, J(mi,m 3 ) - I(m 2 ,m 3 ] 
m\ — m2 



(30) 



I(m 1 ,m 2 ,m 3 ,m 4 ) 



1 


1 


2 2 

mf — m 2 


2 2 

.mf — m 3 



/(mi, m 4 ) 



J(m 3 ,m 4 )^ - 



„ , I(m 2 , m 4 ) - J(m 3 , m A ] 
m 2 — m 3 



(31) 



I(m 1 ,m 2 ,m 3 ,m A ,m 5 ) 
1 



2 2 

mf — m 2 



J (mi, m 3 , 77i 4 , 777 5 ) - J(m 2 , m 3 , m 4 , m 5 ) 



Diagram 1. violation at II and VII from A couplings. 



(32) 



(II, VII) = - £ ^ k A^m ei i mk /(m ei ,m^ k ,m Zm ) + (z <- j 



/,k,m 



(33) 



which can be rewritten as 

(II,VII) = J2^ k ^ kl rn ei rn ek (A + ^teinP)I(m ei ,m^,m ik ) + (i <- j) (34) 
i,k 

when we take A mk to be diagonal. 

Diagram 2. R p - violation at II and VII from A' couplings. 

(II, VII) = -3 £ \' Uk \'^ l m dl A^I(m dl , m^, m Q J + (i <-> j) (35) 

i,k,m 

which can be rewritten as 

(//, V/J) = 3Y / ^ lk ^ M rn dl m d M + ^nP)I(m dn m Dk ,mQ k ) + (i <- j) (36) 

Z,fc 

when we take A™ fc to be diagonal. 

Diagram 3. violation at (IV, VI) in the Grossman-Haber diagram. 

(IV,VI) = a ^ c2 a ( Z **2 ~ Z *ai9' / g) 2 ™ x o a {l(m h , m 9k ,m Dm m x o) cos 2 (a - (3) 

k^-i COS Lj 
, m ' 

+I{m H , m^, m^, m x o ) sin 2 (a - (3) - I(m A , m i>k ,m i)rn ,m x o a 

)} (37) 

Diagram 4. i? p -violation at (I, VII) and (II,VIII) from A' — A' diagram. 

(I,VII) + (II,VIII) = 3 ]T ^^^4^ mdfc A^ mfc if/(m (ifc ,m A ,mQ m ) 

a,l,fc,m Xa 

+3 E ^^^^'^A'^if 

a,l,fc,m Xa 

I(m dk ,m bv mQ m ) (38) 

which simplifies to 

(I,VII) + (II,VIII) = -3j2^^^ ± h k d m 2 dk \' tkk (A + tan/3) 

^K.^^qJ + ^J) (39) 
when both /i^ 1 and A™ 1 are taken to be diagonal. 
Diagram 5 i? p -violation at II and VI from A — A diagram. 

(//, VI) = Y,y j WM jAmlB ^ t ^( 3I ( m eprnE 1 ,m im ,m H+ ) + (i <- j) (40) 

l,m 

which reduces to 



(II, VI) = -J2\ ijl m ej m ei h j e (A + // 4 tan/3)S/tan/37(m e .,m^,m ii ,m ff +) + (i <-> j) (41) 
when we take A mX diagonal. 

1 A 



Diagram 6. 7? p -violation at (I, VII) and (II,VIII) from A — A diagram. 



(I,VII) + (II,VIII) = ]T ^^k h f mek X^A ml I(m ek ,m^m^ 

1 j Tfly/O 

a,l,/c,m Aa 

I{m ek ,mE V miJ 

which simplifies to 



(I,VII) + (II,VIII) = -J2^^^h k e ml k X ikk (A + ^tan(3) 

a,k Xa 

I(m ek ,m Ek ,m ik ) + (i j) 
when both /i^ 1 and A ml are taken to be diagonal. 
Diagram 7. i? p -violation at I and V from A — A diagram. 

(J, V) = E ^-^-^^(^V^ sin 2 P)h^mil(m ej: m H+ ) + (i <-> j) 

a,\,k m X° a 

Taking the lepton Yukawa couplings to be diagonal this simplifies to 



(/, V) = E^^^(^) 2 (^^sin 2 /3)mi/(m ej ,m^ m„ + ) + (i <-> j). 
m xS v2 



Diagram 8. R p - violation at II and V from A — A diagram. 



(//, V) = -E A%^f^ sin 2 (3I(m ej ,m^, m H+ ) + (i <-> j). 



k,m 

which becomes 



(77, V) = -J2 X ijk m ej hi j ^h k e ^ sin 2 (3I(m ej , m^, m H+ ) + (i <-> j). 
for diagonal Yukawa couplings. 

Diagram 9. R p - violation at III and V from A — A diagram. 

(777, V) = E sin 2 ^ nH m ej hi j hfh: k V: 2 U* a2 m+ a 
I(m ej ,m^,m H +,m x +) + (i <-> j) 

which simplifies to 

(777, V) = ]T sin 2 /3 W m e .^ e /i;V Q * 2 C/: 2 m+ 
I(m ej ,mE.,m H +,m x +) + (i <-> ?) 

for diagonal Yukawa couplings. 
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• Diagram 10. _R p -violation at III and VIII from A — A diagram. 



(Ill, VIII) = - £ hfA lk h l rm emf i mH m xt V: 2 U: 2 x 

a,/5,k,l,m 

Zo /34 / ( m e m ,^ k ,m £i ,m+ ) + (i^j) (50) 
when we take diagonal Yukawa couplings this reduces to 



(in, viii) = J2( A + ^ t ^P)KK^m xt v: 2 u: 



a2 >< 



z* z* 

— I(rn ei ,mE i ,m L .,m+ a ) + (i <-> j) (51) 

• Diagram 11. R p - violation at I and VI from A — A diagram. 

(I,VI) = -Y, ^^^KX J rn e] A ml B m t^^I(m e] ,m^m im) m H -) + (i <- j) (52) 

a,l,m Aa 

This reduces to 

(J,\/J) = Yfii a3 a4 (hi j m e .) 2 (A + ^ 4 tan/3)ff 7 tan/3 

I(m ej ,m^,m L] ,m H -) + (i <-> j) (53) 

with diagonal Yukawa couplings. 

• Diagram 12. i? p -violation at III and VII from A — A diagram. 

(///, \//J) = £ fc* i Ik A^ m m em/ i m V;* 2 f/: 2 m x +/(m em , m^ k , m Zl , m,+ ) + (i <-> j) (54) 

a,k,l,m 

With diagonal Yukawa couplings this reduces to 

(///, VII) = - ]T ^m ei (A + m t&nP)X jim m em pL m V* 2 U* 2 m x + 

I ( m e m , m E i , m L i , m x +) + ( i ^ J)' ( 55 ) 

• Diagram 13. _R p -violation at III and VI from A — A diagram. 



(1 1 1, VI) = -^h^A^BitmPh^me^j 

a,k,l 

V a2 U a2 m tJ( m e^ m E^ m L V m H+^ x +) + (i <-> j) (56) 



which reduces to 



Va2K2 m x +H m e 3 ,mE i ,m Li ,m H+ ,m x +) + (i <-> j)- (57) 



1 R 



Diagram 14. i? p -violation at III and IV in g — A loop. 



(III, IV) = -J2 ^KiK2 m x a ^ m ^ h e Bitan/3I{me j ,m Ll ,m H +,m +) + (i <-> j) (58) 

Q,l 

which reduces for g xX — g to 

(JJJ, iV) = - ^ ^U* al V* 2 m+ a ii jmei h ] JBi ttmpi(m ej ,m L ., m H +,m x +) + (i <-> j) (59) 
Diagram 15. i? p -violation at III and VIII in <? — A loop. 

(III,VIII) = -J2 ^ti*m ekl x kN V: 2 U* al m+ a Zhpil^m-^m^) + (i <-> j) (60) 

a/3,k,\ V Z m X0 

When the charged slepton mass eigenstates are aligned with the charged lepton mass eigen- 
states, this becomes 

(///, VIII) = - £ ^X, W lC 2 f/>^/K, ra lo m xt ) + (i <-> j) (61) 
Diagram 16 i? p -violation at I and III in <? — A loop. 

There are two possible diagrams with Iji p at I and III. Firstly the neutralino can arrive as a 
hd at vertex II, where a w~ is absorbed. Secondly the neutralino can arrive as a gaugino, in 
which case the part of the chargino is absorbed. 

g Z* Z* 4 

(I, III) = - -7= sin 2 pWJrrie^iLijV^U^m + pi I(m e . , m H + , m +) 



X /3 



7(m e m H +,m +) + (i <-> j) (62) 



Diagram 17 _R p -violation at I and VII in g — A loop. We also include here the g — A' loop of 
similar form, with sleptons/leptons replaced by squarks and quarks. This diagram, with the 
I/{ p at I removed, corresponds to a mass mixing w°v, or bv. 

We want to show that in the basis where w°v masses are zero at one-loop, the bv mass is finite. 
To do this, we compute both contributions in the = basis. The basis we want is where 
the tree level w°v mass (= g(v i i )/2) cancels the loop w°v mass — so in this basis the bv mass is 
—g'(i>i)/2+ loop = 9 fff 1 ( the w°v loop) + the b°v loop, which is finite. 



w°v mixing gives: 



„k\ y* y* 

(I,VII) = -E^ft^A'^oiCm^miJ 

a,k,l V 2 

-3 £ ^- mdfc ^^^A'^ J B o (0,m (ifc ,m (5i ) + (i <-> j) (63) 
4,1^ m x° a 
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and b — v mixing gives 



(/, VII) 



ak,\ 
ak,l 

+ (*<-► j 







m x° a 1 






\g' kl 


m x° a 


[V2 



'kl 

Qi )-2^X'" kl B (0,m dk ,m 3i ) 



(64) 



Using B (p = 0, mi, m 2 ) = J(m 2 , mi) — \n(m\/Q ) + 1, one can see that the finite contribution 
to b — v mixing, which can contribute to the loop neutrino masses m2 and m\ is the sum of the 
above two expressions (33, 51) with B (p 2 = 0,m 1 ,m 2 ) — > I(m2,m 1 ). 



(Notice that the g' kl appearing in different terms have absorbed different rotation matrices, so 
are not the same. See equation ( |100| ) for formulae with explicit diagonalisation matrices.) 

We can set g kl = g, g' kl = g' and k = 1 in the above expression if the sfermion and fermion 
mass eigenstate bases are the same. 



Diagram 18 R p - violation at III and VII in g — A loop; as noted in ref. [|I2| , this is zero if the 
sleptons are mass-degenerate. 



(Ill, VII) = £ ^m^V^m^Iim^m^m^) + (i <- j) 

which reduces for g kl = g to 

(III, VII) = E^e^ fe K2^ x ±A^/(m efe) m £j) m xi ) + (i <- j) (66) 

a,fe 

Diagram 19 i? p -violation at I and VI in Grossman-Haber loop. 

We neglected this diagram in our previous paper — the expression for the amplitude is lengthly. 
In the interaction eigenstate basis, the incident neutrino can turn into an up-type higgsino/down- 
type higgsino (gaugino), which subsequently turns into a gaugino and up/down Higgs (an up- 
or down-type Higgs and higgsino) at the vertex II. In mass eigenstate basis, these possibilities 
generate a number of terms: 



(65) 



(I,VI) + (IV,VIII) 



1^ n fr — ~m^°(Z 



(32 



Z% 1 g'/g)B l x 



{- \KaL z }-2 - Z^g'/g) sin a + (Z* a2 - Z* a l g'/g)Z} i cos a 

a - f5)I(m h ,m x o g ,m h ] 





+ (z* a2 - 


+ 


ry* ( ry* 

Zj a4\ Zj P2 ' 




+ (z* a2 - 




7* ( 7* 
^4^/32 ' 




+ (z* a2 - 



ay /y^/34 

^i9'/g)coi 
\i9'/g)z^ 
7 4ig'/g)sx 



cos 



+ (67) 
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• Diagram 20 ljl v at I and V with g — A couplings. This differs from figure 2d): the internal line 
is a higgsino and there is an A mass insertion on the scalar line. 

(i,v) = - E ^{(^ + Wff)W + W}^ 

m em hl m fi m sm 2 j3 tan pi(m x +,m^ m ,m H +) + (i <-> j) (68) 
which reduces for diagonal yukawas to 

(i,v) = -T,^§^{(K2 + Ki9yg)u; 2 v; 2 + z: 3 u; 2 v; 1 }m 4 

a, (3 V Xa 

m ej h{nj sin 2 (3 tan (3 x I{rn x + ,m^., m H + ) + (i <-> j) (69) 

• Diagram 21 ^t p at I and IV with g — A couplings. 

The particle identities do not correspond to the labels on figure 2 d). The z/j can mix with a 
higgsino, then the loop fermion is a chargino, arriving at VII as hj,. Alternatively, the z/j turns 
into a gaugino, then there are two possibilities. Firstly, the loop fermion can be a higgsino, 
arriving at II as h u and at VII as ha- Secondly the loop fermion can be a charged lepton, and 
the scalar arriving at VII is a charged Higgs. 



(W) = E ^^{(z: 2 + z: ig '/g)u; 2 v; 2 }m 4 

/3,a,p,m V ^ Xa 

hi m A pm B p t a n(3I{m x +,mz m ,m Lp ,m H+ ) + (i <-> j) 

q,i V 2 ™ x g 
I(m ej ,m~ Lv m H +) + (i <-> j) (70) 



which for diagonal yukawas reduces to 



b? e m e] {A + fa tan (3) Bj tan (3I(m x +,m^., mi.,m H +) + (i <-> j) 

a V2 m xS 
I(m ej ,m L .,m H +) + (i<->j) (71) 



C Basis-independent diagram amplitudes 

In this appendix we write the contributions to the neutrino mass matrix [mjq from each diagram 
in terms of MSSM parameters and invariants. The latter can be evaluated in any basis. 

We list all imaginable diagrams, and explain why some are zero or negligeable. We write the 
MSSM parameters in the mass eigenstate basis, with diagonal quark and charged lepton mass matri- 
ces. See Appendix C for definitions of the matrices Z, U,V,L,E,Q and D which respectively diago- 
nalise neutralinos, negative charginos, positive charginos, doublet charged sleptons, singlet sleptons, 
down-type doublet quarks and singlet down quarks. 



1 Q 



We define 



A ml = -LZ((XA) IM1 ^= + ^X IMI ), (72) 
with L m = X m ■ v/\X m ■ v\, see equation (0). Similarly we define 

Af = -^X'AY ml ^= + ^X' Imi y (73) 

Notice that the indices on A and A d are in the charged lepton and down quark mass eigenstate basis 
respectively. 

C.l A — A diagrams — figure 2a) 

• R p - violation at (I, VIII): this is a loop correction to the tree mass, so can be neglected. 

• R p - violation at (I, VII) and (II, VIII) (diagram 6, equation f42l). 



(J, VII) + (//, VIII) = 5ffi n m e nh™\ii,\A kp E n iE* pl L qm L 

a,l,k,m,p,n,q 



km 



Z '^I(m en ,m E .m i J + (i^j), (74) 



which simplifies to 



(I,VII) + (II,VIII) = -Y8ffi k \v\h k e ml k (A + \Li\Un(3) 



a,k 



Z - ;Z ^/(m efc ,m^,m Zfe ) + (^j), (75) 



rrtyo 

A. a 



when A, L and E are taken to be diagonal. 



i?p-violation at (I, VI) (diagram 11, equation |52|) 



{I, VI) = - Y, S^ihi^A^BME^L^t^P 



Z* 17* 



A.oc 



I(m ej ,Tn^,m Lr ,m H -) + (i^j), (76) 



for diagonal A, L, E this simplifies to 



r7 ^ F~7 ^ 

(I, VI) = Y5i5Uhlm e .) 2 (A+\ii\tan(3)\B\\fi\tan(3 " 3 " 4 

I(m ej} m Sj} m Lj} m H -) + (i j). (77) 

R p - violation at (I,V) (diagram 7, equation |44|) 

(I,V) = -£<5>| 2 (/^mK^m 2 /3W 

a,l,k m Xoc 

+(i^j). (78) 
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Taking A, L and E to be diagonal this simplifies to 

(I,V) = -^^|/,| 2 K^) 2 ^^sin 2 /?tan/5J(m ej ,m,.,m^ + ) + (i <-> j). (79) 



.Rp-violation at (I, IV). This is not possible, because there is no R p - violating mass involving two 
E c . 

violation at (I, II). This is not possible; _R P at II requires three incident leptons. 
i?p- violation at (II, VIII). This is included with R p - violation at (I, VII). 
-Rp-violation at (II, VII) from A couplings (diagram 1 equation 



(II, VII) = - J2 S^S^ n m en A^E np E^L* al L rl I(m en , m Sp , m Li ) + (i «-> j), (80) 

l,m ; n,r,p,q,s 



which can be rewritten as 



(II, VII) =Y,5f k 5{m ei m ek (A + tan/3)/(m ej , m Mk , m Lk ) + (i <-> j), 

l,k 

when we take A, L and E to be diagonal. 



.Rp-violation at (II, VI) (diagram 5, equation ^0|) 



(II, VI) = Sr$B\B\m ej hiA pq tan f3E ns E* qs L* pm L rm 

m,n,s,p,q,r 

I(m ej ,mz s ,m Lm ,m H +) + (i <-> j), (82) 

which reduces to 

(II, VI) = - Sf5 l B \B\m e] m ei hi(A + tan/5) tan/3 

I(m ej ,m Si ,m L[ ,m H +) + (i <-> j), (83) 

when we take A, L and i? diagonal. 
.Rp-violation at (II, V) (diagram 8, equation (46D 



(II,V) =-J2 Sf5^\fi\m e Xrn em tanPsm 2 pE H E* m J(m e] ,mE l ,m H+ ) + (i <- j), (84) 



which becomes 



(/J, 1/) = ~X)^A ^\/j,\m e] hlm e k tan (3 sin 2 (3I(m e] , m^ k , m H +) + (i <-> j), (85) 
for diagonal matrices. 

.Rp-violation at (II, IV) is not possible because there is no _R P mass between two _E C . 
.Rp-violation at (II, III) is not possible; .R p at II requires three incident leptons. 
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-Rp-violation at (III, VIII) (diagram 10, equation |(]) 



(III, VIII) 



a,[3,k,l,m,q,p 



y* y* 

m-yo 



A m e m ^^m iv m xt ) + (i^j), 



when we take diagonal A, L and E this reduces to 



(III,VIII) = J25;5^\ 2 (A+\rit a nf3)(hlml) 2 V: 2 U* a 



a,/3 

y* y* 

Xa 



I(m ei ,m^m- L ,,m x +) + (i <-> j). 



i? p -violation at (III,VII) (diagram 12, equation CT) 



{in, vii) = £ ^r m ^i/^i^^^e m K:2^ 2 ^ x ^n^^,^; fc 

a,k,l,m,n,q,p 

I{rn em ,mE k ,m Li ,m x +) + (i «-> j), 
which becomes, when A, L and i£ are diagonal 



a,m 



i? p -violation at (III, VI) (diagram 13, equation |56|) 



(III, VI) 



which reduces to 



a,k,l,m,q,p 

I(m e .,m^,m l ,m H +,m+) + (i <h> j), 



(III,VI) = ^^NISKA+Mtaa^/i'^me,^ 



I(m e .,mE.,m i .,m H+ ,m x +) + (i j), 



for diagonal A, L and £7. 

-Rp-violation at (III, V) (diagram 9, equation 



(III,V) = £ <^H 2 tan/3r< sin 2 (3m e .hiKE ik E* nk V: 2 U* 



a,k,n 



I(m e rriz ,m H +,m +) + (i j) 
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which simplifies to 



{in,v) = Y, s l s M 2t ^p^ 2 P m e^hiKv: 2 u: 2 m Xo 



I(m e .,m^.,m H +,m x +) + (i^j), (93) 

for diagonal E. 

• Rp- violation at (III, IV), is not possible because there is no ftp mass involving two E c . 

• i? p -violation at (IV , VIII) is not possible, for the same reason as (III, IV) — No diagram with 
E c leaving VII and arriving at V can have i? p -violation at IV, so no diagrams with ft p at 
(IV,...) are possible. 

• violation at (V, VIII) is the same as (I, V). 

• R p - violation at (V, VII) is the same as (II, V). 

• i? p -violation at (V, VI) is not possible, because we cannot put two units of lepton number 
violation on one charged line. 

• .Rp-violation at (VI, VIII) is not possible because H + has no R p conserving mass with E c . 

• i? p -violation at (VI, VII) is not possible because H + has no R p conserving mass with E c . 

.2 A' — A' diagrams — figure 2b) 

• i? p -violation at (I,VIII) is a loop correction to the tree mass. 

• R p - violation at (I,VII) and (II, VIII) from A' couplings (diagram 4, equation |38|). 

(I,VII) + (II,VIII) = 3 £ ?^h k d m d ji: k D kl D* ql Q; m Q rm Af 

a,l,k,m,q,p,r Xa 

I{m dk ,m 6i ,mQ m ) 

+3 £ W^^h^S^DuD^QnJ^ 

a,l,k,m,q,p,r Xn 

I(m dk ,m 6l ,mn), (94) 



which simplifies to 



(I,VII) + (II,VIII) = -S^^H^^^mi^^+Mtan^x 



,k Xo 



I{m dk , mi)k ,mQ k ) + (i <-> j), (95) 
when Ad, Q and D are taken to be diagonal. 

R p - violation at II and VII from A' couplings (diagram 2, equation j35[). 



{II, VII) = -3 E 8l?8$ l m dl A s d r D pk D* rk L nm L* sm 



x 



l,k,m,n,p,r,s 

I{m dl ,m b mn) + {i^ j), (96) 
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which can be rewritten as 



(II,VII) = ^f5^m dl m dk {A+\^nmrn dv m bk ,mQ k ) + (i «- j), (97) 

l,k 

when we take A d , Q and .D to be diagonal. 

Lepton number violation is not possible inside a squark loop, so the only additional effect that 
bilinear Jji v can have is to induce squark flavour violation at V. We neglect this because the 
squark flavour violation is small: 8 l 8f?\n\v u . If we take m u ^ eV then this implies 5^ ^ 10~ 5 , 
so ^li^T^l^u ~ & l x , q GeV 2 , which is negligeable. We differ here from reference [p2 |, where this 



flavour violation due to the bilinears is included, but bilinear lepton number violation is not. 

C.3 the Grossman-Haber diagrams — figure 2c) 

• i? p -violation at (I,VIII) is a loop correction to the tree level mass, so negligeable. 

• i? p -violation at (I, VI) and (IV,VIII) (diagram 19, equation |67]). The basis independent version 
of this diagram can be read off from equation (|57D, substituting 5^\B\ for B n , and 6 l \fj,\ for /i\ 

• i?p-violation at (IV, VI) (diagram 3, equation |37|) 

(IV,VI) = £ 9 ° B ° B J 1 {Z* a2 - Zl l9 ' /gfm x0a LlL ml L) k L nk 

|j(m/ l ,m £ ; ! ,mp fc ,m x c,) cos 2 (a - j3) + I(m H , m vv m Vk , m x o) sin 2 (a - (3) 
-I (m A , m^m^m^ 

)}, (98) 

which for diagonal L becomes 

g 2 5 i <5 J \B\ 2 r 
(IV, VI) = J2 A BB 2 R ( Z *2- z *ai9' / g) 2 m x ° a {l(m h ,m Pi ,m .,m x o a ) cos 2 (a- (3) 

a 4t COS p 

+/(m^,m £;i ,m^.,m x o)sin 2 (a - (3) - I(m A , , m v . , rn x o ) j . (99) 

C.4 g — A diagrams — figure 2d) 

• i? p -violation at (I, VIII) would be a loop correction to the tree level mass, so is negligeable. 

• i? p -violation at (I,VII) (diagram 17, equation |63[). See comments before equation |3. 



(I,VII) = £ ^m efc 5>|^^{25r fc ^^/(^,m e J 

a,k,l,m v Aa 

-25 k x ]m E* kl E m iI(m^ m e j} 
V2 m x o a ■ 



a,kA,m 



-25{* m D* kl D ml I( m£)i ,m dk )} + (i <- j) (100) 

The simplifications when £7, L, Q and -D are diagonal are obvious, so we do not repeat these 
expressions. 
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ftp at (I,V) (diagram 20, equation |B8D . The particle labels in figure 2d) do not apply in this 
case. The incident V{ mixes with a gaugino, and the internal fermion line is a higgsino, or V{ 
mixes with h u and meets w + at II. The scalar incident at VII is an E c . 



m x +m em hism 2 (3 tan (3E j iE* ml I(m x + ,m^, m H +) + (i «-> j). (101) 



This reduces, when E is diagonal to 



%v) = - E «H 2 ^ {(^ + 07#;^ + W^i} «4 

a,p v aq 

m e ,h{ sin 2 /? tan (3I(m + , m©., m#+) + (i <-» j). (102) 



ftp at (I, IV) (diagram 21, equation |7(|). The neutrino can mix with a gaugino, then there are 
two possibilities, which we list separately. Firstly, the internal fermion line can be a higgsino. 
An A insertion on the scalar line ensures that the incident scalar at VII is an E c . The second 
possibility is for the gaugino to interact with £, L at II. In this case H + and e c are incident at 
VII. A third possibility is for the incident neutrino Vi could also mix with h u , which meets a w 
at II, and emits a H£. (An A mass insertion on the scalar line, and m x + on the fermion line, 
ensure that E c and hd arrive at VII. ) This third possibility is a loop correction to the tree 
mass (the loop, without the external mass insertion, is a contribution to /ijZ/j/i°), so we do not 
include it. 



(1,1V), = £ 8^\\B\^= Z «* {Z « 2 + ^ l9 ' /9 \ iE jl E* ml A^L rq L; q tan(3 x 

/3,a,l,m,p,q,r V ^ Xa 

Up 2 V^ m x + I ( m x pmE l ,m Lq ,m H+ ) + (i <-> j). (103) 



The second possibility gives 



(W) 2 = E \-% Z « 3{Z « 2 + Z ^ 9 ' /9) m ej hitzn(3L jm L* nm 

a,m,n V * m Xa 

I(m ej ,m im , m H +) + (i <-> j), (104) 



which reduces, for diagonal A, L and to 



(j,iv) = -Y,tBW\B\^^{(z: 2 + z:j/g)u; 2 v; 2 }x 

m +h J e m e „ (A+ tan/?) tan (3I(m +, m&., mr., m H +) 
+ E WbM I g I -4 K " {Zl2 + Z - MM ™ ej hi tan p 

a m X° a 

I(m ej ,mi m ,m H +) + (i<->j). (105) 



-Rp-violation at I and III (diagram 16, equation |62]) 



a/3 [ m x£ 



-K^^^f^^l + (i - i). (106) 

m x^ 



violation at III and VIII (diagram 15, equation |60|) 



(III,VIII) = - ^^^^L^m^U^m^^^ 

a/3,k,l V ^ '' i X/3 

J(m efe ,m Zi ,m x +) + (z <-> j), (107) 
which reduces for diagonal L to 

(III, VIII) = -Y.^M^KineV^m^ 

i? p -violation at III and VII in g — A loop; this is zero of the sleptons are mass-degenerate 
(diagram 18, equation |65|). 

(III,VII)= £ (^l^^rn^E&m^^ (109) 

a,k,l,m v ^ 

which for diagonal L reduces to 

(/I/, V//) =E*!l/'I^^W^K.'nz il Hi J( j) + (< <- j). (no) 

o,fc 

violation at (III, IV) (diagram 14, equation |58|) 

= -E^4l/i||5|^:i^xi m ^^^tan/? 

I(m ej ,mi l ,m H +,m x +) + (i j), (111) 
which reduces for diagonal L to 

(III, IV) = -^6lS i M\B\^K 1 V: 2 m xt m e .hi^nP 

I(m ep m- Li ,m H +,m x +) + (i <-> j). (112) 

(IV, VIII) and (IV, VII) do not exist, because a doublet slepton must come out of II, so no 
particles carrying lepton number would arrive at VII. 

9fi 



• (IV, V) and (IV, VI) do not exist, because there cannot be two units of lepton number violation 
on a charged line. 

• (V, VIII) and (V, VII) do not exist, because a doublet slepton must come out of II, so no 
particles carrying lepton number would arrive at VII. 



D Conventions and Feynman rules 

We take all coupling constants to be real. Indices on Yukawa couplings and A are in order doublet (d) 
singlet (s): hf s . We work in the mass eigenstate basis of the charged leptons and down-type quarks, 
where the Yukawa couplings are diagonal and only need one index. In Appendix A, we absorb the 
slepton and squark mixing matrices into the Yukawa and gauge couplings at the vertices, so they are 
not diagonal when a scalar meets a fermion. 

The MSSM neutralino mass matrix is diagonalised by the matrix Z m f, with flavour eigenstate 
index /:1..4 corresponding to (—iB,—iW ,h^,h^). The first index m is the mass eigenstate index. 
The chargino mass matrix is diagonalised by matrices U and V: U*MV^ = Md iag , so the positive 
[negative] mass eigenstates are xti = Vmf4>~f [x m — U m fipj] , where ipf = (—iw + ,h+) and ipj = 
(—iw~, hj). The doublet and singlet charged slepton (down-type quark) mass matrices are separately 
diagonalised by matrices Lf m and Ef m (Qf m and Df m ), with index order flavour eigenstate-mass 
eigenstate. We include the A term mixing between doublets and singlets in perturbation theory. The 
matrices L, E, Q and D do not appear in Appendix A; instead we absorb them into A, A^ and the 
gauge and yukawa vertex couplings. We use bold face indices to indicate the squark/slepton mass 
eigenstates in Appendix A, so 

A x ™ = L* pl A™E* qw A l r = Q; i A p fD* qm , h l ™ = L kl h k e m (113) 

We also absorb the diagonalisation matrices into the gauge couplings, so g also appears with indices. 
Usually 

g a = gL* (114) 

however, g appears in the formulae for diagram 17 having absorbed L, E, Q or D. 

We use MSSM Feynman rules from [0], with additional Feynman rules to include the ]ft p interac- 
tions in figures 3,4 and 5. Line direction is superfield chirality. 

For a Lagrangian 

C = $(i$-m)il), (115) 
we fix the phase of the mass insertion to be — % because 

+ t( _ ™)t + - ( 116 ) 



i ~ m if i V 

The — z/ij mass insertions on the external legs of the diagram mix the incident flavour eigenstate 
neutrino z/j with the h° u component of a neutralino. So the external leg propagator delivering a 
flavour eigenstate neutralino / to vertex II is: 

where we have used p 2 = 0. This gives the Feynman rules of figure 3. 

The R p violating /ij and B t mass insertions, and the R p conserving A mass insertion, on inter- 
nal charged lines are negative (by SU(2) antisymmetric contraction), so effectively appear in the 
amplitude with a negative sign: 

; (*N);r A - = 7 1 ~ H ■ (H8) 



•j) — mi "p 1 — rri2 j> — m\^> — rri2 
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(However, the ^ p part of A is — fiiV u / y2 = m ei /Zjtan/?, so appears positive in the amplitude...). 

The Feynman rule we quote in figure 4 for the mass insertion Bi also includes the effect of the 
soft mass m\. The minimisation condition for the potential, in the (Pi) = basis, implies that 
m\ = — Bi tan (3 (at arbitrary loop order) |f[9|| , so the Jft p mixing between H~ and El is 

B- 

- Bi cos 3 + ml sin 3 = *— . (119) 

cosp 



E Numerical Bounds 

We give results for two different values of tan 3 = 2, 10 whenever necessary. This gives us an indication 
of the explicit dependence of the bounds of tan /3, although there is another dependence implicit in 
the neutralino/chargino mixing matrices which we do not address here. Once again, we emphasize 
that we make the assumptions described in section ^ in order that the combination of couplings 
constants are allowed the largest possible values. It is easy to see that under these assumptions the 
integrals that appear in the expressions of appendices A and B can be simply replaced by 



I(m 1 ,m 2 ,m 3 ) 
I(mi,m 2 ,m 3 ,m 4 ) 
J(mi,m 2 ,m 3 ,m 4 , m 5 ) 



167T 2 
1 


777 2 

susy 

l 


167T 2 


m susy 




1 


1 


167T 2 





Below for all bounds we use \m u \. 

• Tree-level contribution 

which gives the constraint 

• Diagram 1 



"V 



Using the mass hierarchy of the lepton sector for i,j^3 we get 



5f 3 S{ 33 < O.leVSTr 2 ^, 

mi 



for m susy = lOOGeV 



(120) 



<l=^>susy, (121) 



6iSi < 10- 12 . (122) 



rrti = -L_L£ 5T k S{ kn m en m ek . (123) 

'"susy kn 



Sf 3 S{ 33 < 2.7 x 10~ 7 , (124) 
For i — 1, 2, j — 3 and j — 1, 2, i = 3 
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Once again using the fermion mass hierarchy we get, 

3h b ml 



Which gives, 



For tan (3 = 2, 
For tan (3= 10, 
Diagram 5 

For tan (3 = 2, 



For tan/3 = 10, 



Diagram 6 



For tan (3 = 2, 



TTi' 1 = 

167i 2 m 



susy 



(^ 3 + ^i? 3 )<0.1evl|^. 

(SIS™ + ^I 33 ) < 3.2 x 10- 7 . 
(« 33 + S'A' 3 ) < 7-2 x 10- 8 . 



lQn 2 m. 



susy 



5f4 





1 


2 


3 


i 





3.2 x 10~ 3 


1.2 x 10~ 5 


2 


3.2 x 10~ 3 





1.2 x 10~ 5 


3 


1.2 x 10~ 5 


1.2 x 10~ 5 






i/j 


1 


2 


3 


1 





1.4 x 10~ 4 


5.5 x 10~ 7 


2 


1.4 x 10' 4 





5.5 x 10~ y 


3 


5.5 x 10~ 7 


5.5 x 10~ 7 






m 



ij — ST 



^ 167T 2 m susy 



kk\ 



sis 



/i A 



i/j 


1 


2 


3 


1 


2.5 x 10~ 5 


2.5 x 10~ 5 


2.5 x 10~ 5 


2 


2.5 x 10~ 5 


2.5 x 10~ 5 


2.5 x 10~ 5 


3 


2.5 x 10~ 5 


2.5 x 10~ 5 


1.0 x 10" 1 



For tan/5 = 10, 



( W + 





1 


2 


3 


i 


5.5 x 10~ e 


5.5 x 10~ e 


5.5 x 10" 2 


2 


5.5 x 10~ e 


5.5 x 10" 6 


5.5 x 10~ 2 


3 


5.5 x 10~ 2 


5.5 x 10~ 2 


2.3 x 10~ 2 



Diagram 7 



,, tan (3 sin (3 „■ „■ . 2 , , 2 2 , , : 
167r 2 m susy " ^ v e > e 



For tan /3 = 2, 



Vj 


1 


2 


3 


i 






7.1 x 1(T 4 


2 






7.1 x 10~ 4 


3 


7.i x lir 4 


7.1 x 10" 4 


3.5 x 1(T 4 



For tan/3 = 10, 



5 l 5 J 

It It 



i/j 


1 


2 


3 


1 






5.7 x 10- 6 


2 




0.35 


5.7 x 10~ 6 


3 


5.7 x 1(T 6 


5.7 x 10~ 6 


2.8 x 10~ 6 



The dashed lines mean that there is no bound. 
Diagram 8 



m j? = V 6 k u m ek Un(3 * in2(3 5f (m e .hi - m e .h[). 



The bounds can be obtained from the bounds of diagram 5/(sin 2 /3). 
Diagram 9 

tan (3 sin 2 /3 



to; 



8n 2 m, 



-5],5{m e .m e .h l p hi. 



susy 



For tan /3 = 2, 





1 


2 


3 


i 








2 






9.3 x 10~ 2 


3 




9.3 x 10~ 2 


3.5 x 10- 4 



For tan/5 = 10, 





1 


2 


3 


i 








2 




0.19 


7.4 x 10~ 4 


3 




7.4 x 10- 4 


2.8 x 10~ 6 



Diagram 10 



The bounds here can be obtained from those diagram 7x tan (3 sin 2 (3. 
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Diagram 11 



to; 



• = tan S i 5 j ( m 2 m 2 ) + SiSUml hf) 

V 167T 2 TO susy f 1 ej e J ^ B\ e t e I 

The bounds here can be obtained from those diagram 7x sin 2 /?. 
• Diagram 12 



(151) 



TO 



k 167T 2 TO susy 



(152) 



The bounds here can be obtained from those diagram 5x tan/3. 

• Diagram 13 

™$ = — ^ 5%{m ei h} e m ej y e ). 

o7T TO susy 

The bounds here can be obtained from those diagram 9x tan (3 sin 2 (3. 

• Diagram 14 

gtan/3 



(153) 



< 3 = f^^W^K + S^ B m ej hi). 



(154) 



For tan (3 = 2, 



W 3 b) 



i/j 


1 


2 


3 


1 




1.2 x 1(T 4 


4.6 x 1(T 7 


2 




6.o x lir 5 


4.6 x 10~ 7 


3 






2.3 x 1(T 7 



(155) 



The % index on 8^ corresponds to the column, and the j index to the row, for both the table 
above and below. 

For tan (3= 10, 



i/j 


1 


2 


3 


1 


0.12 


5.3 x 10~ 6 


2.0 x 10- 8 


2 




2.6 x 10~ 6 


2.0 x 10~ 8 


3 






1.0 x 10~ 8 



(156) 



• Diagram 15 



TO 



'J 



The bounds here can be obtained from those diagram 20 x tan (3 sin 2 /?. 
Diagram 16 

g sin 2 (3 



(157) 



-:^Ji(rn ei h l e + rn e M). 



The bounds here can be obtained from those diagram 20 x tan (3. 



(158) 



^9 



Diagram 17 

From the slepton-lepton loop, neglecting g' with respect to g: 



(5iS{ 3 3 + <^f) = 1.0 x 10- 8 

From the squark-quark loop: 



(^T 3 + ^T) = 3-4 x 1(T 9 

Diagram 18 

This diagram is non-zero only for non-degenerate sleptons. 
Diagram 19 



m 



64-7T 2 cos (3 



in 



susy 



)■ 



For tan (3 = 2, 



W + W = 6.7x10" 

For tan/3 = 10, 



(^b + W = 1-4x10- 

Diagram 20 



,, # tan /3 sin 2 (3 „■ „■ . , , , 



For tan /3 = 2, 





1 


2 


3 


l 




1.5 x 10" 


-4 


5.7 x 10- 


-7 


2 


1.5 x 10~ 4 


7.5 x 10" 


-5 


5.7 x 10- 


-7 


3 


5.7 x 10~ 7 


5.7 x 10" 


-7 


2.8 x 10- 


-7 



For tan/5 = 10, 



if 3 


1 


2 


3 


1 


0.12 


5.4 x 10- e 


2.0 x 10~ 8 


2 


5.4 x 10~ 6 


2.7 x 10~ 6 


2.0 x 10~ 8 


3 


2.0 x 10^ 8 


2.0 x 10~ 8 


1.0 x 10~ 8 
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• Diagram 21 



m » = i^72 {6 ^ Bmejhi + 5 ^ m ^- (169) 

Here we obtain the same bounds as in diagram 14, but with the indices on 5^ and Sb inverted; 
the bound of diagram 14 on S^S B applies to 5^5 l B according to this diagram. 
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Figure 2: Schematic representation of one- loop diagrams contributing to neutrino masses, in a La- 
grangian basis. The blobs indicate possible positions for ]fl p interactions, which can be trilinears (at 
positions II and VII) or mass insertions. The misalignment between jl and v allows a mass insertion 
on the lepton/higgsino lines (at points I, III, or VIII) and at the A-term on the scalar line (position 
V). The soft ]fl p masses appear as mass insertions at positions VI and IV on the scalar line. Figure 
a) is the charged loop with trilinear couplings A (or h e ) at the vertices. Figure b) is the coloured 
loop with trilinear A' or yukawa hi, couplings. Figure c) is the neutral loop with two gauge couplings 
(in the MSSM mass eigenstate basis), and figure d) is the charged loop with one gauge and and a 
Yukawa coupling. This diagram occurs if gauginos mix with charged leptons — that is if 5^ ^ 0. 
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Figure 3: Feynman rules for mass insertions on external legs. The left hand column is a more correct 
representation; we use the abbrieviated notation of the central column in figure 0. 
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Figure 4: Internal charged line mass insertions. 




Figure 5: Feynman rules for trilinear/ Yukawa interactions. For the quarks, replace Xijk — > A' 
E l ~> D%, (e L )j -> (g L ) j; and so on. 
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